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Abstract 

In this paper, we investigate chaotic inflation from scalar field subjected to potential in the 
framework of f(R 2 ,P, Q)-gravity, where we add a correction to Einstein’s gravity based on a 
function of the square of the Ricci scalar R 2 , the contraction of the Ricci tensor P, and the 
contraction of the Riemann tensor Q. The Gauss-Bonnet case is also discussed. We give the 
general formalism of inflation, deriving the slow-roll parameters, the e-folds number, and the 
spectral indexes. Several explicit examples are furnished, namely we will consider the cases 
of massive scalar field and scalar field with quartic potential and some power-law function of 
the curvature invariants under investigation in the gravitational action of the theory. Viable 
inflation according with observations is analyzed. 
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1 Introduction 

The last data mm coming from observations of the universe anisotropy increased the interest for 
inflationary universe. Inflation has been proposed several years ago EHE] to solve the problems 
of the initial conditions of Friedmann universe and eventually to explain some issue related to 
the particle physics, but, despite to the constraints that must be satisfied to fit the cosmological 
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data, the choice of the models is quite large (see Refs. EE] for an introduction to inflationary 
cosmology). 

Many inflationary models are based on scalar field representation, where an homogeneous 
scalar field, called inflaton, is subjected to a potential and produces the accelerated expansion of 
early-time universe, when the curvature is near to the Planck scale. Typically, the magnitude of 
the inflaton is arbitrarily large (chaotic inflation) at the beginning of the inflation |7] and therefore 
the held rolls down towards a potential minimum where acceleration ends: thus, the held starts 
to oscillate and the rehating processes for the particle production take place 0 Ga main], other 
inhationary models are based on a phase transition between two scalar helds Pina- Additionally, 
it is expected that inflation is related with quantum corrections to General Relativity, and in this 
direction many efforts to construct viable models taking into account higher derivative corrections 
to General Relativity emerging at the Planck scale have been done m nsnn ma¬ 
in this work, we would like to investigate how chaotic inflation works in the framework of higher 
derivative corrections to the theory of Einstein. Our modification to the gravitational action of 
General Relativity is expressed in terms of the square of the Ricci scalar and the contractions of 
Ricci and Riemann tensors, in the attempt to include a wide class of models (see Refs. mmm 
for reviews). We will present some explicit examples of theories and potentials for inflaton which 
make viable the inflation according with cosmological data. 

The paper is organized in the following way. In Chapter 2, we present the model: the Hilbert- 
Einstein action of General Realativity is modified by adding a function of the square of the 
Ricci scalar and the contractions of Ricci and Riemann tensors, and the contribute of a scalar 
field subjected to a potential is also included. Thus, we derive the Lagrangian and the related 
Equations of motion of the theory in flat Friedmann-Robertson-Walker space-time, by using a 
method based on Lagrangian multipliers which reduces the Equations of motion at the second 
order. It is interesting to note how for Friedmann-Robertson-Walker metric such a kind of models 
has an equivalent description as Gauss-Bonnet theory. In Chapter 3, we investigate the general 
feauture of chaotic inflation from scalar held in the framework of our modified theory. Chapter 
4 and Chapter 5 are devoted to the explicit examples, corresponding to two well-know cases of 
chaotic inflation, namely chaotic inflation from massive scalar held and chaotic inflation from held 
with quartic potential. In the framework of General Relativity they lead to viable inflation, and 
we are interested to see how results change for some toy model of Gauss-Bonnet modified gravity 
and some model based on the power-law functions of the curvature invariants under investigation. 
The conclusions with some final remarks are given in Chapter 5. 

We use units of kp = c = h = 1 and denote the gravitational constant, Gn, by k 2 = 8itGn, 
such that Gn = 1/Mpj, Mp\ = 1.2 x 10 19 GeV being the Planck mass. 


2 Formalism 


Let us consider the following action 

R 


1= / V~9 
J M 


2r2 + f(R 2 , P, Q) - - V(</>) 


k 2 = 


Stt 

Wi 


( 1 ) 


where A4 is the space-time manifold, g is the determinant of the metric tensor g^ v and R is the 
Ricci scalar. The gravitational part of the Lagrangian takes into account the higher derivative 
corrections to Einstein’s gravity encoded in the generic function f(R 2 , P,Q) = f where 

P = R^RT , Q = R^R^ , (2) 


Rfj, v and R^a^ being the Ricci tensor and the Riemann tensor, respectively. The “matter” part of 
the Lagrangian depends on a scalar field <f> subjected to the potential 

We will work with the flat Friedmann-Robertson-Walker (FRW) metric, whose general expres¬ 
sion is given by 

ds 2 = — N(t) 2 dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ), y/—g = N(t)a(t) 3 , (3) 
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where TV(t) = TV is a lapse function and a(t) = a is the scale factor, the both depending on the 
cosmological time t. Thus, the curvature invariants of the model on flat FRW space-time read 


R =-^A X + Y )’ P =^( X2 + Y 2 + XY ) > Q=^( x2 + Y2 )’ ( 4 ) 


TV 2 


TV 2 


with 


X----— V- — 

a aN ’ a 2 ’ 


( 5 ) 


and the dot denotes the derivative with respect to the time. By plugging this expressions into the 
action m, we obtain an higher derivative Lagrangian. However, by using a method based on the 
Lagrangian multipleyer EH eh [H2 [22, we can deal with a first order standard Lagrangian. We 
introduce the Lagrangian multipliers £ , a , £ as 


I = 


Na* 




k 2 N 2 


(X + Y) + f{R\P,Q) - C ( P- {X + Y) 


12 


12 


-°[ p -n*( x2 + y2 + xy ))-z[Q-m-A x2 + Y2 ))+^ 2 ~ V W 


TV 4 


2TV 2 


( 6 ) 


In order to get a first order Lagrangian, we rewrite the action as 


I = 


Na* 


IM 


k 2 N 2 


(X + Y) + f(R 2 ,P,Q)~((R~ — (X + Y) 


+ 20 + ^-v'W , (T) 

such that the variations with respect to R , P and Q lead to 

c = /p(P 2 , P,Q) + =f- [.fp(R 2 , P, Q) + /q(P 2 , p, Q)] , <7 = /p(P 2 , P,Q), £ = /q(P 2 , P, 0), 

( 8 ) 

where /p,p,q(P 2 , P, Q) are the derivatives of /(P 2 ,P, Q) respect to R,P and Q, 

,2 D rf/(P 2 ,P,Q) , , d2 D ^_ df{R 2 , p ,Q) f , d2 Dn , d/(H 2 ,P,0) 


/«(P 2 ,P,Q) = 2P- 


, / P (P 2 ,P,Q) = 


, /q(P, P, Q) = 


dP 2 ’ 

Thus, after integration by part, we obtain for the gravitational part 
(a,a,N,R,R,P,P,Q,Q) = - ^ + /TVa 3 - CPiVa 3 - 


dQ 


( 9 ) 


—o' P — 


R“ 

P- 

and the total Lagrangian results to be 

fc'tot — fc'grav H - £<j) i ^(f) — 


TVa 3 - £ Q - 


R' 


Na 3 + -jp (d + 2^ , (10) 


2TV 


<^ 2 - V(<f>)N a' 


( 11 ) 


As a result, we obtained a first order Lagrangian respect to the unknown variables TV(t), a(t ), P(f) = 
R,P(t) = P,Q(t) = Q. 

Some remarks are in order. Obviously, the expression m can be generalized to the case 
f(R 2 ,P,Q) —> f(R,P,Q). An interesting special case is given by Gauss-Bonnet gravity. The 
Gauss-Bonnet four dimensional topological invariant reads 

24 

G = P 2 - 4P + Q, G = —XY , (12) 
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where the second expression is the form of the Gauss-Bonnet on the flat FRW space-time. If 
f(R, P, Q) = f(R, G), by taking into account that 

MR, p, Q) -f MR, G)+2RMR, G ), MR, p, Q) ->■ -4 MR, G ), f Q (R , p, q) -t MR, G ), 

(13) 

we derive from OB, 

iV, R, R, 0,6) = -2^ + JV„»(/ - W f ° • < 14 > 

according with Ref. [23]. Moreover, it is possible to demonstrate that the Lagrangian of f(R, P, Q)- 
models corresponds to the Lagrangian of f(R,G)~ theories on FRW background [251126] , We may 
replace Q with Q = G — R 2 + 4 P and therefore / q(R , P, Q) with /g(R, G, Q ) in (HOD and make 
the following substitutions 

MR, P, Q) fn(R, G, Q) + 2Rf G (R, G, Q ), MR, P, Q) fp(R, G, Q)-Af a (R, G, Q ), (15) 

in order to cancel the additional derivatives that we have aquired. We get 

£ grav (a,a,N,R,R,G,G,Q,Q) = — ^° 3 (/ ~ Rfn ~ Gfa) — 

-fp( P + Na 3 M + 


N 


N 


and after integration by part we obtain 


Co 


,(a,a,N,R,R,G,G,Q,Q) = + Na 3 (f - Rf R - Gfa) ~ 


+fpNa 3 


n 2 N 


R 2 


6 aa 2 f R 

6 da 2 /it 

8d 3 ■ 
iV 3 * G 

N 

N 

+*(£ 

4 a 2 aR' 

) ,(16) 

N 

6 aa 2 f R 

6 aa 2 f R 

8d 3 ■ 
iV 3 ^ G 

N 

N 


G 


3 ~ ~2 P 


(17) 


but on FRW metric the last term is null and we recover (Hill . To pass from f(R 2 ,P, Q) to f(R, G)- 
gravity on FRW space-time, we can substitute R 2 ,P ,Q with R 2 ,G ,C 2 into the action 0. Here, 
C 2 is the “square” of the Weyl tensor, 


/~<2 _ 1 p2 o r? PM*' i D 

O - 3 -Tl ZlLf^iyll + Hi, ±\ , 


and the following relations are met, 


R 2 


M2 M E> 2 

P "-2 + T’ 


(18) 


(19) 


On FRW metric ([3]), the square of the Weyl tensor is identically null ( C 2 = 0, 6C 2 = 0) and does 
not contribute to the dynamics of the model, such that we can drop down it from the Lagrangian 
and use the formalism of f(R 2 , G)-gravity. 

By making the variation of m cu respect to N(t) and therefore by putting N(t) = 1, we 
find 

^ + (/ - c R)+6H 2 (+6H(-v (p - y) -Z (q - m)~ uh3 + 2 ^) = Y + V{ ^ ’ (20) 

where H = a/a is the Hubble parameter. The variation respect to a(t) with N(t) = 1 leads to 

±(3H 2 + 2h) +(/-C-R) + (47?C + 6iL 2 c) + 4H( + 2( - a (p - ^ - £ (q - 

/ o \ 

(a + 2 $J 


-8H (H 2 + H} (& + 2^j 


— AH 


= - y-^) 


( 21 ) 
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Finally, the variations respect to R,P and Q for the gauge N(t) = 1 read 

R = 12H 2 + 6 H, P = 12 [h 2 + 3 H 4 + 3 HH 2 } , Q = 12 (h 2 + 2H 4 + 2 HH 2 ^} . (22) 

In conclusion, we obtained a set of five second order differential equations EDI) (122[) . By taking 
the time derivative of EDI) and therefore by using ETTl we derive the energy conservation law of 
the field, 

4> + 3 Hj> = -V'(</>). (23) 

Here, the prime denotes the derivative with respect to <fr. Thus, we can make the following 
identification, 

P<!> = Y + v {4>) > P<t> = y - v {4 >), H + 3i?(^ +P0), (24) 

where and p$ are the effective energy density and pressure of the field, respectively. We also 
may introduce the Equation of State (EoS) parameter of the field as 


u = p ± = 2 V (</0 

- P<t> <i> + 2V{<P) ' 

Let us see now how the inflationary cosmology is reproduced by such a kind of models. 


(25) 


3 Inflationary cosmology 

We would like to see how a modification to Einstein’s gravity based on the higher derivative 
correction terms R 2 ,P ,Q changes the classical picture of the inflation from scalar field models in 
the General Relativity background. The dynamics of the model © is governed by the equations 
(EDI) and EDI) with (EDI) . 

The inflation is described by a quasi-de Sitter solution, with the Hubble parameter which 
slowly decreases with the time, such that the slow-roll approximations are valid, 


, H , 

I I 1 ’ 


1-^, 

HH 


« 1. 


(26) 


It means, that the magnitude of the so called “slow-roll parameters” 


H 

IP ’ 


rj = 


H H 
lp ~ 2 HH 


^ = 2e- 


1 


2eH ’ 


(27) 


must be small during inflation. Moreover, 0 < e in order to have H < 0, and, since the acceleration 
is expressed as 

- = H + H 2 , (28) 

a 

we see that the accelerated expansion finishes only when the e slow-roll parameter is on the order 
of the unit. 

To describe the early-time acceleration, we will use the approach of chaotic inflation. At the 
beginning the field, namely the inflaton, is assumed to be negative and very large. The slow roll 
regime takes place if the kinetic energy of the field is much smaller respect to the potential, 

0 2 <C V(0), \4>\ <C 3 Hj> (29) 

In this way, the field EoS parameter E5l) is ~ — 1 and the de Sitter solution can be realized. 
On the other hand, from EDI) we have, in the slow-roll approximation EDI) . 

3 3 H4>~-V" {(/>)<]), (30) 
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and, if V'(a) > 0, the kinetic energy increases with the field which tends to a minimum of the 
potential, where EH is not valid and inflation ends. 

For the (quasi) de Sitter solution of inflation, one introduces the e-folds number as 


M = log 


Of 

Oi 


H(t)dt ~ 3 


r<t> i 


H 2 


' <t> f 




dtp, 


(31) 


where a\j are the scale factor at the beginning and at the end of inflation, the related times and 
<pi t { the values of the field at the beginning and at the end of inflation. Generally, the primordial 
acceleration can solve the horizon and velocities problems of Friedmann universe if 55 < A f. 

By using the slow-roll parameters, one can also evaluate the universe anisotropy coming from 
inflation deriving the spectral indexes. The amplitude of the primordial scalar power spectrum 
reads 


A 


2 _ 
■r ~ 


k 2 H 2 

&7T 2 € ’ 


(32) 


and according with cosmological observations must be ~ 10 9 ; the spectral index n s and the 
tensor-to-scalar ratio r are given by 


n s = 1 — 6e + 2rj , r = 16e . 


(33) 


The last results observed by the Planck satellite |2] constrain these quantities as n s = 0.9603 ± 
0.0073 (68% CL) and r < 0.11 (95% CL). 

For chaotic inflation in the background of General Relativity, namely in the case of action 0 
with f(R 2 , P , Q) = 0, one has in terms of the field potential and its derivative, 


i (v'M V 

2k 2 V V(<P) ) 


d= — 

K z 




where the quasi de Sitter solution of inflation is given by H^ s = n 2 V(tp)/ 3, while, in the slow-roll 
approximation, H ~ k 2 V'(tp)<p/(6H<is) with cp derived from (1501) . 

In our case, the quasi de Sitter solution of inflation H — H^s, H^s being a constant, is given 
by equation (T201) under the condition (EH1) . namely 


+ (f(R 2 s ,P ds ,Q d s ) 

R'dsf p (RdS’ -^dS; Qds) _ 

4 


Rdsfli(RdSi PdS,Qds) ^ 


V(tP), 


RdsfQffidS’ RdS ? Q ds) 

6 

(35) 


with 

R dS = 12H 2 s , P dS = ^ , Qds = ^- (36) 

By taking the derivative of (EH with (1221) and by using the slow-roll conditions (EH and EH . we 
obtain the equations for H , H : 


H 


6 H 


+ l2Hf R - lMH 3 f RR - 5184 H r f PP - 2304 H 7 f QQ 

— V(tp)<P > 


1728 H 5 f RP - 1152 H 5 f RQ - 6912 H 7 f PQ 

(37) 


H 


6 H 


+ 12Hf R - 144 H 3 f RR - 5184 H 7 f PP 


2304 H 7 f QQ - 1728 H 5 f RP 


~2V\cp)'<p-2HV\(p)Ppe, e 


H 

IP’ 


1152 H 5 f RQ - 6912 H 7 f PQ 


6 


(38) 

















where H = H^s and </> is determined by (1551) . The last term of (1551) has been approximated to 
simplify it: it comes from terms proportional to ~ H 2 , which are negligible if potential is flat 
(in such a case, e <C |r/|), but not in other cases like for the power-law potentials that we will 
analyze. Thus, the dependences on V((/)), V'(cf>) and V"(<j>) of the slow-roll parameters e , 77 and of 
the e-folds number change respect to inflation in Einstein’s background case. 

In the following chapters, we will consider some suitable potentials for scalar field and we will 
see how the early-time acceleration is reproduced according with cosmological data in /(f? 2 , P, Q)- 
gravity. 


4 Quadratic potential 


One classical example of chaotic inflation is given by massive field with potential 

V{4>) = 0 cm, (39) 

where m is a positive mass term smaller than the Planck mass, m -C Mpi, in order to avoid 
quantistic effects during inflation. Let us assume <f> negative and very large. In the slow-roll 
regime (1551) . equation (1501) leads to 


log 


3PdS 


(t — ti) + const, 


(40) 


where Pas is the quasi-de Sitter solution of inflation given by (1551) and </>; the value of the field at 
the beginning of inflation when t = tp For example, in Einstein’s gravity with /(P 2 ,P, Q) = 0, 
one has 


PdS 



/ 7 r m(f >i 
3 Mp\ 


3« 2 {t tl) ~ 


rriMpi 

v'm 


(t - U). 


Thus, the slow-roll approximation (1291) is valid as soon as 


(41) 


Mpi 

\/l27r 


< \<k\ < 



(42) 


where we have also taken into account that H^s < Mpf. however, the accelerated expansion ends 
only when the e slow-roll parameter is equal to one. Note that the field is larger than the Planck 
Mass, but its kinetic energy is smaller. The e-folds and the slow-roll parameters (1541) read 


M* Pl Mp t 2^ 2 

47 uj)? ’ 47 T(j)? ’ Mpi 


(43) 


and for large e-folds (i.e. (f >i much larger than the Planck Mass) the slow-roll parameters are small 
and the spectral index n s in (1551) can satisfy the Planck data. However, we must stress that the 
tensor-to-scalar ratio r results to be slightly bigger than the Planck bound. In general, this is true 
for all the power-law potential in scalar field representation (except for the case V (</>) ~ <j> , which 
is negative for large values of the field and presents some criticism). In the last months, the correct 
value of the tensor-to-scalar ratio has been a debated question, and for this reason the analysis of 
such a kind of models is still interesting. For example, the last BICEP2 results mi indicated for 
the P-mode polarization of the CMB-radiation the tensor-to-scalar ratio r = 0 . 201 q'q 5 ( 68 % CL), 
and the vanishing of r has been rejected at 7.0er level, while, as we stated before, the data coming 
from Planck experiment reveal an upper bound for the tensor-to-scalar ratio at r = 0.11 with 95% 
CL. Moreover, it should important to mention that very recent combinations of the Planck and 
revised BICEP2/Keck Array likelihoods lead to r < 0.09 with 95% [28] . 

Let us see now how inflation induced by massive scalar field works in some toy models of 
/ (i? 2 , P, Q)-gravity. First of all, we will consider the subclass of Gauss-Bonnet gravity, and then 
we will extend our investigation to more general theories. 
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4.1 Gauss-Bonnet models 

The Gauss-Bonnet (1331) in the flat FRW space-time ([3]) with the gauge N(t) = 1 is given by 


G = 24 H‘ 


(H 2 + H) . 


For the case f{R 2 , P, Q) = /(G), by taking into account (1351) . equation (1331) reads 


3 H 2 


+ (/ - G/g) + 24 H 3 Gf GG = r- + V(<t>) 


( 44 ) 


(45) 


k‘ " 2 

where we have used the fact fa = G/gg- Thus, the de Sitter solution is the corresponding of (1351) . 


3 H; 


dS 


+ (/(G dS ) - G dS /G(G dS )) = V (</>), G dS = 24 H\ s , 

rx 

and by taking the derivative of (1351) one has in the slow-roll approximation (1351) , 

6 HH 


- 2304 HH 7 f GG ~ V\(t>)<t >, 

!X~ 

which corresponds to (1371) with (1T31) . Moreover, the equation for H is derived as 


6 HH 


- 2304 HH 7 f GG - 2V'(<j>)4> - 2 e , e = - 


H 

IP' 


Let us introduce the quadratic potential (1331) and assume the following form for /(G), 

/(G) = jG n , n^l, 


(46) 


(47) 


(48) 


(49) 


where 7 is a dimensional constant such that [ 7 ] = [ Mp^ and n is a number. For n = 1 we 

recover Einstein’s gravity being the Gauss-Bonnet a topological invariant in four dimension. 

The simplest non trivial case of (1331) is given by n = 1/2, for which the dimension of 7 is [ 7 ] = [M|> ; ]. 
We may assume 0 < 7 and introduce an effective mass of the theory as 


TO e ff = m. 


3M pi 


3 Mp t + 87t-\/67 


, 0 < 7 , TOeff < TO , 


such that, in analogy with (1311) (1331) . one has the following solution of (1351) . 

pK TOeff^i . . , TOeffMp; , 


H dS = -2 


3 Mpi ’ 


vT 2 tt 


-{t - U ), 


and 


Mpl < |^| < M2pl 


~ 0.5- 


. M 2 


pi 


Vl2TT ' '' rl ' ' 2TOeff V 7T TO e ff 

The slow-roll parameters (1371) and the e-folds (15T1) follow from (1371) (l48l) as 

2 , ,.9 / \ 2 


(50) 

(51) 

(52) 


e ~ 


Mk(PlL) Ml (IM ) 2 /cox 

47T0? \ TOeff / ’ 47 T(j)f \ TOeff / ’ Mp; V TO / 


and for m e ff = m we recover (1351) . Thus, for large boundary values of the field </>;, the e-folds Af 
can be large enough and the slow-roll parameters are very small during inflation, since 
















like in the case of /(G) = 0. 

The amplitude of the primordial scalar power spectrum usa of the model is 

2 _ 16 mj s 7 T (/> 4 

n 3 M^m 2 ' 

The spectral index and the tensor-to-scalar ratio (EH read 


M 2 Pl . 

( m ) 2 


( m\ 

tt;/? 

\m e sj 

n(/)f 

\m e gj 


and in order to satisfy the Planck data (see under Eq. (1551) 1 one must require 

0.0324 < f) < 0.0470, 

7r</>f \m eS J 


( 55 ) 


(56) 


(57) 


or, in terms of the e-folds number, 

42 < Af < 61. (58) 

In the range 55 < J\f < 61 inflation is considered viable (we also must stress that acceleration 
continues after the slow-roll approximation and ends only when e = 1). Thus, 


2.6 

m eS 


< \<t>i\ <3.1 


mMpi 


m e ff 


(59) 


Despite to the fact that the spectral index n s is viable, the model presents the same criticism of 
the massive scalar field in General Relativity framework concerning the tensor-to-scalar ratio r, 
which results to be r ~ 0.13, slightly bigger than the Planck constrain at r < 0.11. We note that, 
since the curvature during inflation cannot exceed the Planck mass, the condition (l59l) with (1521) 
leads to 

m < 0 . 2 Mpi , (60) 

and in order to recover the primordial scalar power spectrum (l55l) cs 10 9 , it must be 
to ~ 10 ~ 6 Mpp The Gauss-Bonnet contribution /(G) = 7 \[G to the action is compatible with the 
inflationary scenario from scalar field with quadratic potential, but leads to a tensor-to-scalar ratio 
slightly bigger than the one given by the Planck data. If 0 < 7 , the inflation is realized at curvature 
smaller respect to the classical case with 7 = 0 (on the other side, if —3M|, i /(87rv / 6) < 7 < 0, we 
obtain the opposite behaviour), but a suitable setting of the initial value of the field permits to 
recover the same spectral index. 


Let us take now 1 < n in (l49l) : in this case, the value of /(G) ~ 7 R 2n is dominant respect 
to the Hilbert Einstein term during inflation when 

(iff) 2 " _1 <R< M2pl ’ (61) 

such that the equations (1551) (1551) can be solved by neglecting the contribution coming from R/n 2 . 
The de Sitter solution exists and it is real if 7 < 0, 

HdS - $(-</>) 1/2n , /’^/’i + (-^i)^^r(i-£i), $ = (— = m , (62) 

3$ y v / 2(24) n 7(l - n)) 


with [$] = [Mpi] . The slow-roll parameters (1571) and the e-folds number (1511) read 


6?r4> 2 (—</>i) - 


V - 


3$ 2 (—0i)l 


a r~ 


3n(—</>,) ™ T 2 


(63) 
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For large values of fa the e-folds is large and the slow-roll parameters small. However, we will see 
that the Planck constraints limit the magnitude of fa. We observe 


1 n 

6 ~ W v ~ Jr ’ 


such that e < 77 . The amplitude of the primordial scalar power spectrum (1321) is given by 

6n(—</>i) 2 / n <b 4 


A 2 - 


mPMpi'K 


and the spectral index and the tensor-to-scalar ratio pID are derived as 

(3 — 2 n)m 2 8 m 2 

Tl ^ - - rp — - 

3n4> 2 (— fa)i ’ 3?r4> 2 (— fa)i 

In order to satisfy the Planck data we must require 

0.0324 < ( 3 ~ 2n ) w < o 0470 l < n < |. 
3n$ 2 {-fa)i 2 


(64) 


(65) 


( 66 ) 


(67) 


The condition on n is quite restrictive and imply that inflation described by the model takes place 
very near to the Planck scale. For example, if n ~ 3/2“, it follows from (RTT1) . 


-7= < R < Mpi , ( 68 ) 

but in this case the magnitude of the boundary value of the field in (1671) cannot be very large and 
the e-folds of the model is quite small. A better fit of the cosmological data may be found for a 
value of n between one and 3/2, but the inflationary scenario produced by the model cannot be 
defined “chaotic” due to the restrictions on the boundary of the field. 

Viable inflation based on the account of qG" , 7 < 0,1 < n, could be realized only if 1 < n < 
3/2, but, as soon as n is close to 3/2, the magnitude of the field is small even if the curvature is 
extremely near to the Planck scale. On the other hand, if n is close to one, condition (RJT1) is not 
well satisfied and the model turns out to be the one with /(G) = 0. In the last part of the next 
subsection, we will reconsider such a model by adding a contribution from R n : we will see that 
also in this case the conditions on n for a feasible inflation do not change. 


4.2 f(R 2 ,P,Q) power-law models 

In this subsection, we will consider an explicit model of f(R 2 ,P,Q) in the context of chaotic 
inflation from quadratic potential. To simplify the problem, we will rewrite P, Q as funcions of 
the square of the Ricci scalar P 2 , the Gauss-Bonnet invariant G and the square of the Weyl tensor 
G 2 as in (fTTHl . Since the contribution of the Weyl tensor is identically null on FRW metric, we can 
reduce the theory to f(R 2 , G)-gravity, and the equation (f^TJl) with (ITT?1) reads 

^r + if-Rfu- Gfa) + 6 H 2 fn + 6 Hf R + 24 H 3 f G = ^ + V(fa). (69) 

Thus, the de Sitter solution is derived from 

+ (/(/&> Gas) - R d sfR(R 2 dS ,G dS ) - G d sf G (R 2 dS ,G dS )) + 6 H 2 s f R (R 2 s ,G dS ) = V{fa ), 

AC 

P dS = 12U 2 S , G dS = 2AH\ S . (70) 

The derivative of (1M1) in the slow-roll approximation (l^()l) . 

^- + 12 HHf R + UAH 3 Hf RR - 1152H 5 Hf RG - 2304 HH 7 f GG ~ V’(</>)$ , (71) 

Av 
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6 HH 


corresponds to (1371) with (THU) . and the equation for H is given by 


+12 HHf R + U4H 3 Hf RR - 1152H 5 Hf RG - 2304 HH 7 f GG - 2 2 HV'(<j>)j> e, e = - . 

(72) 

For our purpose, let us take the following Ansatz for f(R 2 ,P,Q)~ model, 

f(R 2 1 P,Q)=aR 2n + pP m +7Q p , (73) 

where ck,/3, 7 are dimensional constant and n,m,p numbers. We get from (1191) . 

f(R\G 1 C 2 )=aR 2n +p(^-^ + ?pj + 7(2 C 2 -G+^pj . (74) 

The behaviour of this model on FRW space-time correponds to the behaviour of 

*- + ,(*-f)\ T (*- C )’, (75) 


being the Weyl tensor identically zero on FRW metric. To study how inflation can be realized 
from such a kind of theory, we must do some assumption. For m = p = 1, the model reduces to 
f(R 2 )- gravity, being G a topological invariant in four dimension (its contribution in (RTH) drops 
down), and we get 

f(R 2 ) = aR 2n +tR 2 , m = p = 1. (76) 

For n < 1, we deal in fact with a R 2 correction to standard gravity. In this case, at subplanckian 
scale, the Hilbert Einstein term R/k 2 is dominant in the action, and inflation has a corresponding 
(viable) description in the so-called Einstein frame [2Hj after a conformal transformation of the 
metric. In literature we have many studies about inflation from R 2 in Einstein frame [30] . or 
inflation from R 2 combined with other curvature invariants coming from trace-anomaly, quantum 
corrections or string inspired theories mm- For 1 < n we obtain more general power-low cor¬ 
rections to General Relativity m ■■ note that if we neglect the Einstein’s term we do not have real 
de Sitter solution for positive values of the potential, and n must remain close to one. 


One simple non trivial case is given by n = m = p = 2 in (1731) (1751) , for which we get 


f{R 2 , G) = aR 4 + j3G 2 + jR 2 G , £ = (a - 


n = m = p = 2 , 


where 


, /3 7 n P , ~ 

° = a+ 9 + 9- /, = I+T’ 7 = -3~T' 


with [a] = [/3] = [ 7 ] = [1/MpJ. Inflation takes place in high curvature limit, 


M 2 Pl 


< R < M 


pi 1 


(77) 

(78) 

(79) 


where <5 is a term with the dimension and magnitude of a , (3 , 7 : in this case the Hilbert-Einstein 
contribution can be neglected in the action. The de Sitter solution is derived from (1701) as 


Has = 4>(—</>i) 1/4 , ^~<^i + (-0 i )4^-(t-t i ), $ = 


i ) 1 / 4 


3$ 


( 2 7 / 83 i/ 4 ) (-365 - /3 - 67 ) 


1/8 


(80) 


with [<f>] = [Mpj 4 ] and (365 + /3 + 67 ) < 0. The behaviour of the field, as usually, is governed by 
m- The slow roll parameters and the e-folds number read 


m 2 (365 + ft + 67 ) 
1203 + 37 - 725)\/^$ 2 




3V^$ 2 ’ 


M ~ 




m 


(81) 


11 









In order to have e > 0 with a real solution for the de Sitter solution (1551) . it must be 


P + %1 _ (P + 67) 

72 a 36 

From m we get 

_ (36a + p + 67) _ _2_ 

6 ~ 20 + 37 - 72a)U ’ V ~ N ' 

The amplitude of the primordial scalar power spectrum (1321) is derived as 

A 2 _ 320 + 37 — 72d)(—</>j)<l> 4 
U m 2 Mp l 'ir( 36d + P + 67) 

and for the spectral index and the tensor-to-scalar ratio (1551) one finds 


(82) 


(83) 


(84) 


n. ~ 1 — 


m 2 (396d — P + 67 ) 


60 + 37 - 72 < 5) v /^$ 2 

To satisfy the Planck data it must be required 

0.0324 < ra 2 (396d — /3 + 67 ) 


4m 2 (36d + P + 67 ) 
r = - 


3(/3 + 37- 72 


6(/3 + 37- 72d)V 3 W$ 2 


< 0.0470, 


(85) 


( 86 ) 


but also in this case the tensor-to-scalar ratio r is bigger than the Planck contraints, being on 
the order r ~ 0.26. In order to recover ~ 10~ 9 , it is enough to have m ~ 10 ~ 6 Mpi. We 
immediatly see from (1551) . 


a < 


P - 67 

396 


(87) 


Conditions (1551) and (1571) must be satisfied simultaneously. We have several possibilities. If 7 = 0, 
namely /3 = — 2 y in (1751) (1751) . 


^ <0 - 4 < “ < 556 ’ ’ = °- < 88) 

but in this case the theory is affected by antigravitational effects during inflation, since the effec¬ 
tive gravitational constant of the model, G e s = Gjv/(1 + 2 K 2 fn(R 2 ,G)), Gn being the Newton 
constant, results to be negative if a < 0. 

In general, if 7 < 0 and 67 < ft < — 67 , conditions (1551) and (1571) can be satisfied for positive 
values of a with the possibility to avoid antigravitational effects. By using A f of (1811) in (1551) we 
get 

21 0 + 37 — 725) <u< 31(/3 + 3-y — 72a) 

(396d — ft + 67 ) (396d — P + 67 ) ’ 1 

and to have 55 < W we must require (P + 37 — 72d) ~ 2(396d — p + 67 ), such that a suitable 
value of pi can reproduce a sufficient amount of inflation according with Planck results. 

We have seen that the model f(R 2 , P, Q) = aR A + PP 2 + 7 Q 2 with a massive scalar field may 
bring to a viable early-time acceleration with curvature near to the Planck scale, but the tensor- 
to-scalar ratio is not compatible with the Planck data (but, for example, may be compatible with 
the BICEP2 results [27]). When the curvature decreases, we are out of the range (1791) and the 
Hilber-Einsten term becomes dominant in the action: at this point, the reahting processes for 
particle production take places and Friedmann expansion starts. 


To conclude this chapter, we would like to reconsider the model (l49l) with 1 < n of the previous 
subsection togeter with a power law function of the Ricci scalar in the context of /(i? 2 , G)-gravity, 
namely 

f(R 2 , G) = aR 2n + 7 G" , 1 < n, (90) 
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where a, 7 are constants whose dimensions are [a] = [ 7 ] = [ 1/Mp[ n 1 ' ) ]. Inflation starts at high 
curvature, 

<r<k '' <9i) 

where <5 is a term with the dimension and the magnitude of a , 7 , such that the Hilbert Einstein 
contribution R/k 2 can be ignored into the action. In the presence of massive scalar field, the de 
Sitter solution of the model reads 


(92) 

with [$] = [Mp™ 1 ^" n ] and ( 6 ™a + 7 ) <0. Therefore, the slow roll parameters and the e-folds 
number are derived as 

m 2 (n — l)(6™a + 7) 


6n((n — 1)7 — 6 n na)(— ^i) 1 /”$ 2 ^ 3(—</>i) 1 /” < I> 2 

In order to get e > 0 with a real solution for the de Sitter solution we must find 


( 93 ) 


(n - 1)7 a _J_ 

6 n n 6" ‘ 

We also observe 

(n — l)( 6 n a + 7 ) n 

2((n — 1)7 — 6 n na)J\f ’ ^ A f 

Thus, the amplitude of the primordial scalar power spectrum (1321) is derived as 

2 6 n( 7 (n — 1 ) — 6 n na)(— 0 i ) 2 / n $ 4 

n m 2 Mp l n(n — l)( 6 n a + 7 ) 


(94) 


(95) 


(96) 


and for the spectral index and the tensor-to-scalar ratio we get 

, ?u 2 (6 n (n(3 + 2n) — 3)a + (n(5 — 2n) — 3 ) 7 ) 8m 2 (n — l)( 6 "a + 7 ) 

^ ^ — _______ j* — _._ 

3n((n — 1)7 — 6 n na)(—(f>i) 1 / n ^> 2 ’ 3n((n — 1)7 — 6 n na)(— ()>i) 1 /™<I > 2 

(97) 

As a consequence, to reproduce the Planck data we must require 


0.0324 < 


m 2 (6"(n(3 + 2 n) — 3)a + (n(5 — 2 n) — 3)7) 
3 n((n — 1)7 — 6 re ua)(-((ii) 1 /"$ 2 


< 0.0470, 


(98) 


and finally 


(3 — (5 — 2n)n)"f 
6™(n(3 + 2 n) — 3) 


(99) 


In order to satisfy conditions (11741) and (1M1) . we can take 7 < 0 and 0 < a (avoiding antigravitational 
effects) only if 


3 

1<H< 2 ’ 


( 100 ) 


recovering the same result of (1(771) where the R 2n contribution was not considered. In other words, 
the addition of a R 2n contribution to the model in (737)1) does not change the range of n, which 
remains quite close to one to reproduce the Planck results. 
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5 Quartic potential 

Let us consider now a quartic field potential in the general action 

y ($ = ^r> 0<A ’ ( 101 ) 

where A is a positive adimensional constant. In Einstein’s gravity where f(R 2 1 P, Q) = 0, for large 
and negative value of the field we get from (1551) and (15U1) . 


#dS 


[2tt VX $ 

V T Mpi ’ 


Vam 


PI 




(t - ti) . 


( 102 ) 


where, as usually, </>i < 0 is the boundary value of inflation. Thus, the slow-roll approximation 
(1551) is valid as soon as 

tB < l<Ai| < Mpi Ox) ’ (103) 

where Afp;/(A ) 1 / 4 -C 1 and the field may be larger than the Planck Mass during inflation. In this 
case, the e-folds and the slow-roll parameters (1551) read 


e ~ 


M P , 


V ~ 


3Mp, 

2t T<f>? 


Af ~ 


M 2 Pl 


(104) 


and, for large e-folds, the slow-roll parameters are small and the spectral index n s in (1331) can 
satisfy the Planck data, with a larger value of the tensor-to-scalar ratio r. 

In the following subsections, we will see some significative examples of f(R 2 ,P,Q )~gravity 
where chaotic inflation from field with quartic potential could be realized. 


5.1 Gauss-Bonnet models 

As a first example, we will consider the Gauss-Bonnet model 

/(G) = 1 VG, 0 < 7 , (105) 


7 being a positive dimensional constant such that [ 7 ] = As we have already seen in S 14.11 

this kind of correction to Einstein’s gravity leads to a viable inflation in the presence of massive 
scalar field. 

If we introduce the effective A e ff parameter, 


A e ff = 


A Mp t 


M 2 pl + (4tt/3)7\/24 


, 0 < 7 , Aeff A , 


for the potential (1 101 1) . Equation (1551) leads, on the de Sitter solution, 

^27T \J Aeff 4> 2 ill VA esMpi 


with 


H dS = 

Mpi 

y/Sn 


3 Mpi 


< \<j>i\ < Mpi 


\Z6tt 


(t - U) , 


3 A 1/4 _ 1.5 Mpi 

27T Aeff' / 


A. 


1/4 ’ 

eff 


(106) 


(107) 


(108) 


in analogy with (11021) (11031) . 

The slow-roll parameters (1571) and the e-folds (1551) are derived from (1571) (1551) and result to be 


~ 7T(j) 2 y A e ff' J ’ 


(±) 

27 T(/>? \XcbJ 




7T0f ( Aeff A 

M 2 Pl \x)' 


(109) 
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and for A e ff = A we recover HMD. Since 


3 


( 110 ) 


1 

Jr' 


we see that for large boundary values of the field </>,, the e-folds A f can be large enough and the 
slow-roll parameters very small during inflation. 

The amplitude of the primordial scalar power spectrum (13211 of the model is given by 


A 2 27r A~ff0f 

n 3M%\ ’ 

while the spectral index and the tensor-to-scalar ratio (1551) read 

•- wtf \ >,,.»] ■ ntf 1 ,'VI. / ■ 

Thus, in order to satisfy the Planck data we must find 


( 111 ) 


( 112 ) 


0.0324 < < 0.0470, 

7 T<pf \Aeff / 

or, in terms of the e-folds number, 

64 < AT < 93 . 

As a consequence, the boundary value of the field must be 


4.5M P; 


't— <\4>i\ < 5.4 Mpi 

Aeff 



(113) 


(114) 


(115) 


As in the case of scalar field with quartic potential in the framework of General Relativity, the 
tensor-to-scalar ratio r is bigger than the Planck bound (r ~ 0.16). Since the curvature during 
inflation cannot exceed the Planck mass, the condition above with (I108D leads to 


Vts <0 - 33 ' (I16) 

and in order to obtain ~ 10 -9 , one must have ^/A/VC ~ 10~ 7 . In conclusion, the model 
/(G) = 7 \/G ,0 < 7 , in the presence of scalar field with quadratic (see § 14.11) or quartic potential 
may lead to a viable inflationary scenario, but the tensor-to-scalar ratio r exceed the Planck result. 
Moreover, in the case investigated in this subsection the e-folds in (I114D is larger than the e-folds 
in (1551) . As a consequence, the slow-roll parameters of the model are smaller in the presence of 
scalar field with quartic potential respect to the quadratic potential case. 

5.2 f(R 2 ,G) power-law models 

In this last subsection we will discuss a general form of /(i? 2 ,G) power-law model and we will 
see that, in the presence of scalar field with quartic potential, differently to the cases of quadratic 
potential analyzed in § 14.21 inflation cannot be realized. We also remember that, as we have 
already discussed, a /(i? 2 , G)-model can be seen as a representation of a f(R 2 , P, Q)-theory on 
FRW space-time if we use (flT)l) and therefore consider null the Weyl tensor and its contributes to 
the action. 

Let us start by the folliwing f(R 2 , G)-model, 

f(R 2 ,G) = aR 2n +l3G n +j(R 2 G) n/2 , 1 <n, (117) 
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a,/3 ,7 being dimensional constants such that [a] = [f3] = [ 7 ] = [Mp , 1 n) ]. This correction to 
Einstein’s gravity is dominant respect to the Hilbert-Einstein term R/n 2 in the action at high 
curvature when 

(if ) 15=1 < R < M ™ ■ (118) 

where S is a term with the dimension and the magnitude of a , / 3 , 7 . The de Sitter solution of the 
model for the quartic potential (I101D follows from (1701) as 


H dS = <S>(-fa) 1/n , $ = 


^1/(4 n) 


2(3n+2)/(4n)3l/4( n _ ]_)l/(4n) (_6 "q. _ fj — 6™/2^ 


1/ (4n 


(119) 


with [<I»] = Mp™ 1 ^ n and 6 "a + /3 + 6"/ 2 7 < 0. Thus, the slow-roll parameters (1071) and the e-folds 
m are derived as 

2(n — 1) ( 6 ”a + /3 + 6 n / 2 7 ) A (—_ x{-fa^-^ 3 $ 2 

3n((2/3 + 6 re / 2 7 )(n-1)-2"+ 1 3 ra ?za) $ 2 ’ V ~ $ 2 ’ ~ ^(n-i) ’j ' 

(120) 

We immediatly see that, for 1 < n, large value of fa leads to large slow-roll parameters and small 
e-folds, rendering the inflationary scenario unrealistic. 

This result cannot be considered a general behaviour of f(R 2 , P, Q)- or f(R 2 , G)-gravitational 
models: however, we would like to note that our Ansatz (1TT71) repersents a quite generic and 
reasonable power-law model of f(R 2 , G). For such a form of correction to Einstein gravity we can 
say that chaotic inflation from scalar field does not work in the presence of quartic potential. 


6 Conclusions 

In this paper, we have investigated chaotic inflation with scalar held subjected to potential in the 
framework of f(R 2 , P, Q)-modihed gravity, namely the gravitational action of the theory includes a 
correction based on an (arbitrary) function of the square of the Ricci scalar R 2 and the contractions 
of Ricci (P) and Riemann ( Q ) tensors. This form of modified gravity is quite general, and the 
curvature invariants under consideration may be related with quantum corrections to General 
Relativity or string inspired theories. To derive the Equations of motion on flat Friedmann- 
Robertson-Walker space-time we used a method based on Lagrangian multipliers and we treated 
the curvature invariants as independent functions: as a consequence, we deal with a system of 
second order differential equations simplifying the analysis of the model, which leads in principle 
to fourth order differential equations. We note that on FRW metric every f(R 2 ,P, Q)-theory can 
be reduced to a Gauss-Bonnet f{R 2 ,G)- theory, since in fact one of the curvature invariant can 
be expressed in terms of the other two. This feauture is manifest by replacing the contractions of 
Ricci and Riemann tensors P, Q with the Gauss-Bonnet G and the square of the Weyl tensor C 2 : 
on FRW metric the Weyl tensor and its derivatives disappear and we can drop down its contribute 
from the Lagrangian. We used the f(R 2 , G)-representation to study our models, since in this way 
the Equations of motion result to be simplified. 

Chaotic inflation from scalar field with potential can be realized in the framework of higher 
derivative models as well as in the framework of General Relativity, but, despite to the fact 
that the continuity equation of the field keeps the same form in the two theories, the Hubble 
parameter and its derivatives depend on the field potential in different ways. Inflation must 
satisfy several constrains to be “viable”: the (quasi) de Sitter solution must take place and the 
slow-roll approximations must be valid, it means, the slow-roll parameters must be small, the e- 
folds sufficiently large to guarantee the thermalization of observed universe, and the spectral index 
and the tensor-to-scalar ratio must satisfy the Planck data. We presented the general formalism to 
investigate inflation and its characteristic paramters in f(R 2 , P, Q)-gravity with scalar held and we 
furnished several explicit examples. In the specific, we investigated two well-known forms of chaotic 
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inflation, namely chaotic inflation from massive scalar field (quadratic potential) and chaotic 
inflation from field with quartic potential. We confronted the results in the framework of General 
Relativity and in the one of our modified theory giving some examples of corrections to Einstein’s 
gravity based on power-law functions of the Gauss- Bonnet or of the other curvature invariants 
under investigation. The (positive) corrections based on the square root of the Gauss-Bonnet are 
on the same order of magnitude of the Hilbert-Einstein term in the action at high curvature and 
permit to realize an inflationary scenario similar to the Einstein’s case. More interseting are the 
higher power-law functions of the Gauss-Bonnet and the other curvature invariants. In this cases, 
at high curvature the modification to gravity is dominant with respect to the Hilbert-Einstein 
term in the action and drives inflation. Thus, by fitting the parameters and the boundary value 
of the field, we may recover a viable inflation in the case of scalar field with quadratic potential, 
but with quartic potential the inflationary scenario appears to be unrealistic. We stress that this 
result cannot be read as a general behaviour of /(i? 2 ,P, Q)- or /(R 2 , G)-gravitational theories, 
even if our Ansatz for the presented power-law models is quite generic and reasonable. Moreover, 
even in the presence of our Ansatz, we cannot state that this models are not able to reproduce 
an early -time acceleration in agreement with Planck data, but only that in the context of chaotic 
inflation induced by large magnitude value of the inflaton inflation is not viable. 

A last remark is in order about the tensor-to-scalar ratio number of the models investigated 
in the present work, which results to be larger than the Planck contrain. This feauture is quite 
general in chaotic inflation from power-law potential, but, since the exact value of the tensor-to- 
scalar ratio is still object of a debated question, we think that this kind of models have to be still 
investigated. On the other hand, the attempt of our study is to furnish a general formalism for 
chaotic inflation in higher derivative gravity theories, which is valid independently on the specific 
examples here analyzed. 

Other works on higher derivative corrections to Einstein’s gravity, FRW /(G)-gravity and 
inflation can be found in Refs. [5511541155] . 
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